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Abstract : We study the normalized eigenvalue counting measure da of 
matrices of long-range percolation model. These are (2n+l) x (2n+l) random 
real symmetric matrices H = {H(i,j)}ij whose elements are independent 
random variables taking zero value with probability 1 — if) {{i — j)/b), b G M + , 
where if) is an even positive function if)(t) < 1 vanishing at infinity. It is 
shown that if the third moment of ybH{i,j), i < j is uniformly bounded 
then the measure da := da n ^ weakly converges in probability in the limit 
n,b — > oo, b = o{n) to the semicircle (or Wigner) distribution. The proof 
uses the resolvent technique combined with the cumulant expansions method. 
We show that the normalized trace of resolvent g n ,b(z) converges in average 
and that the variance of g n ,b{z) vanishes. In the second part of the paper, 
we estimate the rate of decreasing of the variance of g nt b(z), under further 
conditions on the moments of ybH(i,j), i < j. 

AMS Subject Classifications : 15A52, 45B85, 60F99. 

Key Words : random matrices, semicircle law, percolation model. 
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1 Introduction 



Spectral theory of random matrices is a relatively new branch of mathe- 
matics that intensively develops due to its rich mathematical content and 
aho due to numerous applications in theoretical phisics, wireless communi- 
cations, financial mathematics and other fields (see reviews monographs [5], 
[T3] . [T9] . [20] and references therein). 
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In theoretical physics random matrices of infinitely increasing dimensions 
are used to give the statistical description of systems with large number of 
degrees of freedom. 

The first application was in nuclear physics, where E.Wigner proposed 
to consider the eigenvalues of real symmetric random matrices as a model 
for the energy levels of heavy atomic nuclei (see the monograph [IH]). The 
real symmetric random matrix A n of size 2n + 1 introduced by E. Wigner is 
defined by 

A n (i,j) =iV- 1 / 2 a(*,j), \i\,\j\ <n, (1.1) 

where N = 2n + 1 and {a(i,j), i < j} are independent and identically 
distributed random variables defined on the same probability space (Q, P) 
such that 

E{a(i,j)} = 0, E{a 2 (t,j)}:=v 2 , (1.2) 



where E is the mathematical expectation with respect to P. 

Denoting by A_ n < ... < An the eigenvalues of A n , the normalized 
eigenvalue counting function is defined by 

a n {\,A n ) := N-^{\f ] < A}. (1.3) 

E.Wigner [2T] proved that in the case when a(i,j) has all moments finite, 
the eigenvalue counting measure da n (\,A n ) weakly converges in average as 
n — > oo to a distribution g?<t sc (A), where the nondecreasing function cr sc (A) is 
differentiable and its derivative p sc writes as follows 

f V4t; 2 -A 2 -r |\| < 9 7 , 

MA)=.„(A) = { - ,f ll-^ e . 

This limiting distribution (1.4) is referred to as the Wigner distribution, 
or the semicircle law. 

Since then, the convergence to the semicircle distribution was proved 
for various random matrix ensembles that generalize the Wigner ensemble 
[H [12l [16]. Among them we cite the band and the dilute random matrix 
ensembles M, M, CDS EH] . 

In the band random matrices model, the matrix elements take zero value 
outside of the band of width b n along the principal diagonal, for some positive 
sequence of real numbers. This ensemble can be obtained from A n (1.1) by 
multiplying each a(i,j) by I(-i/ 2 ,i/ 2 ) ((i-j)/b n ), where 



lA(t) 



1 if t e A, 
if t e R \ A 
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is the indicator function of the interval A. The ensemble of dilute random 
matrices can be obtained from A n (1.1) by multiplying a(i,j) by independent 
Bernoulli random variables 



d{i,j) 



1 with probability p n /n, 

with probability 1 — p n /n, < p n < n. 



Assuming that b(n) = o{n) for large n, the semicircle law is observed, after 
proper renormalization,for both ensembles, in the limit b n — > oo [H] and 
p n — > oo [9] as n — > oo. 

It should be noted that the matrix {d(i,j)}ij represents the adjacency 
matrix of random Erdos-Renyi graphs (see [6] for definitions and more de- 
tails). Random Erdos-Renyi graphs are extensively studied. Recent applica- 
tions also involve the problems of spectral theory of random graphs models 
[7] and references therein. 

Our work is motivated by papers where the random graphs of certain 
percolation models are considered. More precisely, we are interested in a ge- 
neralization of the tow ensembles mentioned above. Roughly speaking, we 
consider band random matrices with a random width. To proceed, we mul- 
tiply each matrix element a(i,j) by b n d(i,j), where < b n < N is a 
sequence of integers and d(i,j), \j\ < n is a Bernoulli random variable 
with 

1 with probability i>{—' 



(o with probability 1 ^' 5) 

and < if)(t) < 1 is an even function vanishing as t — > oo. We assume that 
T^n = {d(i,j), \j\ < n} is a family of independent random variables, 

also independent from A n = {a(i,j), \j\ < n}. 

The family {d(i,j), \j\ < n} can be regarded as the adjacency 

matrix of the family of random graphs {T n } with N = 2n + 1 vertices 
such that the average number of edges attached to one vertex is b n . Hence, 
each edge e(i,j) of the graph is present with probability ip((i — j)/b) and 
not present with probability 1 — ip ((i — j)/b). Below are some well known 
examples : 

— In theoretical physics, the ensemble {T n } with ip(t) = e~'*' s is referred 
to as the Long-Range Percolation Model (see for exemple [T] and re- 
ferences therein). Our ensemble can be regarded as a modification of 
the adjacency matrices of {T n }. To our best knowledge, the spectral 
properties of this model has not been studied yet. 

— It is easy to see that if one takes b n = N and = 1, then one recovers 
the Wigner ensemble (1.1). 



3 



— If one considers (1.5) with ip(t) = I(-i/2,i/2)(t), one gets the band ran- 
dom matrix ensemble described before [T5] . 

The paper is organized as follows. In section 2, we define the random ma- 
trix ensemble H n of long-range percolation model and state our main results. 
In section 3, we use the cumulant expansions method to derive relations for 
E{g n (z)}, where g n (z) is the Stieltjes transform of the normalized eigenvalue 
counting measure of our random matrix ensemble. In section 4, we show that 
Var{g n (z)} vanishes and prove the semicircle law for the random matrix H n 
under the condition that supj • E{|o(i, j)| 3 } < oo. In section 5, assuming more 
conditions on the moments of a(i,j), we derive more precise and powerful 
estimates for the variance Var {g n ( z )} ■ 

2 The ensemble and main results 

Let us consider a family of independent real random variables A n = 
{a(i,j) : |j| < n} satisfying (1.2) and defined on the same probability 
space (0,F,P). 

Let ip(t), t e R, be a real continuous even function such that : 

< ij)(t) < 1 and / if>(t)dt = 1. (2.1) 
Jr 

Introduce a family of independent Bernoulli random variables V n = 
{d(i,j) : |j| < n} as in (1.5), that is independent of the family A n . 
We assume that A n and V n are defined on the same probability space. 

Define a real symmetric N x N random matrix H n by : 

H n (t,3)=b- 1/2 a(t,3)d(i,j), i<j, \i\,\j\<n, (2.2) 

where < b n < N, N — 2n + 1. We study asymptotic spectral properties of 
the ensemble {H n } in the limit when n — > oo with 

lim b n = oo and lim — = 0. (2-3) 

n— yoo n— >oo n 

This limit (2.3) when 1 << b n « n is known as the limit of relativity narrow 
band width. 

Here and below the family {H n } is referred to as the ensemble of random 
matrices of long-range percolation model. 
Let us introduce 

/2 r = supE{|a(z,j)r}, r e N, (2.4) 

hi 
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the absolute moments of a(i,j). 



Our main result is as follows. 

Theorem 2.1 // fa < oo (2.4), the normalized eigenvalue counting 
function a n (X,H n ) in (1.3) for the ensemble given in (2.2) with ip satisfies 
(2.1), converges in probability in the limit (2.3) to a nonrandom distribution 
function : 

p- lim a n (\,H n ) = a sc (X), (2.5) 

n— *oo 

where cr sc (X) is given by (1.4). 

Let {/x n } n be a sequence of probability measures that converge weakly to 
a probability measure fi. Let g n and g denote the Stieltjes transforms of \x n 
and \i respectively defined by 

f fj. n {d\) f fadX) 
9n{z) = / ; g(z) = / , (2.6) 



A — z J A — z 

for z G C \ R. Then the weak convergence of /i n is equivalent to the conver- 
gence of the Stieltjes transform (in the uniform topology on compact sets 

Let cr n (dX, H n ) be counting measure of the eigenvalue of H n : 

a n (d\,H n ) = ^Y, 6 ^ dX ^ ( 2 - 7 ) 

|i|<n 

where 8\ t is the Dirac mass at Aj. Then Theorem 2.1 is equivalent to the 
following statement. 

Theorem 2.2 Let g n {z) be the Stieltjes transform of the normalized ei- 
genvalue counting measure a n (dX, H n ) (1.3). If fa < oo (2.4), then in the 
limit (2.3), there exists a nonrandom function w(z) such that for all nonreal 

z 

p - lining g n (z) = w sc (z), (2.8) 
where w sc (z) is the Stieltjes transform of the semicircle law given by (2.5). 

Remarks 

1. The limiting function w sc (z) is the solution of equation 



that maps C + — > C + when 

C+ = {z e C : Imz > 0} 

and the parameter v is determined by (1.2). 
2. To prove (2.8), we start by proving it for z G A,,, where 

A, = G C : 7] < \Im(z)\} with rj = 2v + 1. (2.10) 

The convergence in probability (2.8) is shown in a standard way : 
first we prove that E{g n (z)} converges to w sc (z) (2.9) and that *Var{g n (z)} 
vanisfes. More precisely, we are going to prove the following statement. 

Theorem 2.3 If fa < oo (2.4) and z E A v (2.10), then in the limit 
n — > oo (2.3), 

lim B{g n (z)} =w sc (z) (2.11) 

n— >oo 

cmc? 

Var^)} = Oib- 1 ' 2 ), (2.12) 
where w sc (z) is given by equation (2.9). 

Theorem 2.2 follows then using Tchebychev's Inequality and gives Theo- 
rem 2.1. 

To prove Theorem 2.3, we mainly follow papers pUI IE]- We apply the 
cumulant expansions method of [8] and get similar relations to the ones 
obtained in [10]. So, combining these two approaches, we develop a method 
to study the resolvent of H n . It is rather general and can be used to study 
other random matrix ensemble with jointly independent entries. 

It is not hard to show that the estimate (2.11) together with (2.12) imply 
convergence in probability (2.8). We give the details at the end of section 3. 
Then (2.8) implies convergence (2.5) ( see section 3 for the proof). 

Then in section 3, we give the limiting equation from E{g n (z)}, In section 
4, we prove Theorem 2.3 and we deduce from this the results of Theorem 2.2 
and Theorem 2.1. In section 5, we give more details on asymptotic behavior 
of the variance and prove auxiliary statements. 

For the sake of simplicity we omit the subscript n in b n , G n and H n , so 
that G n = G (2.12) is the resolvent of the matrix H n = H (2.2). 
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3 Limiting equation for ~E{g n (z)} 

The resolvent 

G n {z) = (H n - zl)~\ Imz 0, 

is widely exploited in the spectral theory of operators. Its normalized trace 
coincides with the Stieltjes transform g n (z) (2.6) of the normalized eigenvalue 
counting measure <r n (d\,H n ) (1.3) ; 



[ cr n (d\,H n ) J_ >^ 1 

9n[Z) ~ J X-z ~ N X {n) 



\j\<n A j Z 

Let us consider the resolvent identity for two hermitian matrices H and 
H' : 

G-G' = G(H - H')G', (3.1) 

where G{z) = (H - zl)' 1 , G'(z) = {H' - ziy 1 . 

Regarding (3.1) with H = H n (2.2) and H' = 0, we obtain relation 

G(i,j) = £8v-ZY< G &P) H b,j), Z = -z-\ (3-2) 

\p\<n 

where 5 denotes the Kronecker symbol 



1 if i — j 
if i^j, 

H(i,j) = H n (i,j), \j\ < n are the entries of the matrix H n and G(i,j) 
are the entries of the resolvent G n = (H n — zl)^ 1 . 
Relation (3.2) implies that 

EG(M) = E{G(i,p)H(p,i)}. (3.3) 

\p\<n 

To compute the mathematical expectation E{G(i, p)H (p, i)} , we use the cu- 
mulants expansion method proposed in [8]. 



3.1 Cumulant expansions and principal equation 

Let us start with the description of our basic technical tools. 

(z). The cumulant expansions formula. Let us consider a family {Xj, j = 
1, .., m} of independent real random variables defined on the same probability 
space such that E{|X,| 9+2 } < oo for same q 6 N and j = 1, ..,m. Then for 
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any complex-valued function F(ti, .., t m ) of the class C(R m ) and for all j, one 
has 

Efoi^, .,X m )] = ± ^n dTF{X d l - Xm) } + , (3-4) 

r=0 J 

where K r is the r-th cumulant of X, [8j and the remainder e q can be estimated 
by the inequality 



|^|<Csup |^^|E{|X J f+ 2 } (3A) 

and C is a constant depending only on q. 

Relations (3.4) and (3.5) can be proved by using the Taylor's formula (see 
e.g. [8] for example). At the end of section 5, we prove (3.4) in the case of 
q = 5 and give more details on the form of the remainder terms e J q of (3.4). 

It is known that the cumulants can be expressed in terms of the moments 
of Xj. If we denote E(X j ) = E(Xf) = E(X|) = and fx r = E(XJ) with 
j = 1, ..,m, one obtains : K\ = [i\ = 0, K 2 — H2, K3 = 0, X 4 = /i 4 — 3^, 
K 5 = 0, Kq = [Mq — 15/^4/^2 + 30/x|, etc. 

(ii). If if is a real symmetric NxN matrix with elements H(i,j), < 
n, and G = (H — z)~ l {i, j) is its resolvent, then 



dG st _ \ -G sj G kt si j = k 



dH jk \ -G sj G kt - G s kGjt si j ^ k. 
The relation (3.6) follows from the resolvent identity (3.1). 



(3.6) 



E{G(i,p)H(p, i)} = -V{^g^M^) + e ip . 



Now, let us explain the main idea of the proof of relation (2.15). Applying 
(3.4) to the resolvent identity (3.3), one gets 

Substituting this equality in (3.3) and using (3.6), we arrive at relation 
EG(i, l ) = Z + ^-J2 V{G(t, i)G(p, p)M 1 -^) + R n ,b{i), 

\p\<n 

where R n ,b(i) vanishes as n, b — > 00 (2.3) (see subsection 3.2 for more details). 

If one assumes that the average E{G(i, z)6 _1 ^ G{p,p)ip(\p — £]/&)} fac- 
torizes, then g n ,b{i) = EG(i,i) satisfies equality 

g n ,b(i) = £ + —g n ,b(i) ^2 9n,b(p)il>(—g-) + RnA*)- 

\p\<n 
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Assuming that the value b" 1 ^2 p ip([i — p\/b)g n ,b(p) does not depend on i in 
the limit n,b — > oo, we get convergence g n ,b(i) —> w(z), for all \i\ < n — bL, 
where L is sufficiently large and w(z) is the solution of equation (c.f.(2.9)) 

w{z) = i + iv 2 w{ z )\ i = - z -\ 

3.2 Derivation of relations for EG(i, i) 

Let us consider the average E{G(i, p)H (p, i)} . For each pair (i,p), G(i,p) 
is a smooth function of H(p,i). Its derivatives are bounded because of equa- 
tion (3.6) and the inequality 

I^,P)|<||G||<^, (3.7) 

which holds for the resolvent of any real symmetric matrix. Here and below 
we denote by ||e||| = ^M<n l e WI 2 the euclidean norm of the vector {e(i)}|j| 



i\<n 



and \ \G\ \ = sup|| e || 2=1 HGe]^. 

Inequality (3.7) implies that, \D pi Gi P \ — M\Imz\~ 3 where M is an abso- 
lute constant. In what follows, we use the notation D pi = d r /dH r (p,i). 

According to (2.2), (1.5) and the condition fi 3 < oo, the third absolute 
moment of H(p,i) is of order b~ 3 ^ 2 . Applying (3.4) to E{G(i,p)H(p,i)} with 
q — 1, we get relation 



E{G(i,p)H(p,i)} = K&{G{i,p)) + K 2 E(DlMhp)) + e ip (3.8) 
where K r is the r-th cumulant of H(p,i) and satisfies the estimate 

\e ip \ < ^su P \D 2 pi G(hp)m\a(p,i)\ 3 }E{\d(p,i)\ 3 } 

-W^^ E{I ° (P ' ?)|3}E{I ^' ?)|3} ' (3 ' 9) 

for some constants C and M. 

To compute the partial derivatives, we use (3.7). For a complex symmetric 
matrix G, one has D£G(i,i) = (-l) a a\G(i,i) a+1 , with a = 1,2 and 

Dlfi^p) = -G(t,p) 2 -G(t,t)G(p,p), (3.10) 
D 2 pi G(i,p) = 2G(i,pf + 6G(i,t)G(p,p)G(i,p), (3.11) 
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for distinct i and p. Using formulas (3.4) and (3.5), the cumulants can be 
expressed in terms of the moments. Considering X = H(p,i), we get 

h = ^ = E(X) = ±=E{a(p,i)}E{d(p,i)} = 0, (3.12) 

K 2 = & = E(X 2 ) = jE{d(p, if} = (3-13) 
Using (3.10), (3.12) and (3.13), we rewrite (3.8) in the form 

2 

E{G(i,p)H(p,i)} = - V -E{G{i,pf + G(i,i)G(p,p)}i>(^-) + e tp , p + i 

(3.14a) 

E{G(i, i)H(i, %)} = - V -E{G(i, t) 2 M0) + e u . (3.146) 

o 

Substituting (3.14) in (3.3), we obtain equality 

EG(M) = * + ^ £ E{G(.,p) 2 M^) 

\p\<n, p^i 



|p|<n, pj^i \p\<n 



Relations (3.9) and (3.11) imply the following estimate 



(3.15) 



|c *' ^ |7^v^ ( — } (3 - 16) 

for some constant C. 

For a given a random variable, we write (° for : 

c° = c-ec, 

so that (3.15) takes the form 

EG(M) = £ + yEG(M) £ EG(p,p)^(^) 

|p|<n 

|p|<n 
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where 

<f>i(i) = E E{G(*,p) 2 M^) - ^E{G(M) 2 M0), |i| < n, 

\p\<n 

(3.18) 

^« = X^ E{G(? ' ?)G ° (p ' p)Mi T Z) ' (3 ' 19) 

\p\<n 

where G° = {G(i,j) — E(G(i, j))}ij and e ip is given by (3.16). We are inter- 
ested in the average value of the normalized trace of the resolvent 

g n (z) = N- l Y J G(i,i), 

\i\<n 

where G(i,i), \i\ < n are the diagonal entries of the resolvent. Taking the 
normalized sum of (3.15) and using the notation g n (z) = g n ,b(z), we obtain 

E VG(i,i)VG(p,p)^-> 

\i\,\p\<n 

+^E^«}+^E^)}-| E ( 3 - 2 °) 

|i|<n |i|<n KUp|< n 

To prove limE{(? ni 6(z)} = (2.11), we need the following statement 

concerning the pointwise convergence in average of the diagonal entries G(i, i) = 
G nt b(i,i; z), \i\ < n of the resolvent. Given a positive integer L, let 

B L = B L (n,b) = {i e Z : \i\<n-bL}. (3.21) 



Lemma 3.1 Given e > 0, there exists a positive integer L = L(e) such 
that 

sup \E{G n , b (i,i;z)}-w(z)\ < e, z e A v (3.22) 

ieB L 

for sufficiently large b, n satisfying (2.3). 

Lemma 3.1 will be proved in the next subsection. Let us assume that 
(3.22) is true. Then regarding definition of Bl, we deduce from inequality 
(3.7) that 

1 1 — n+bL— 1 1 n 

ie{^)}-^E eg mi = i^ E eg m + ^ E eg mi 

ieB L i=~n i=n-bL+l 
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2bL-2 , . . 

< T, 2 G A„. 3.23 

Now (2.11) follows from Lemma 3.1, inequality (3.23) and the limiting condi- 
tion (2.3). o 

3.3 Proof of Lemma 3.1 

Let us have a step back and rewrite (3.17) as 

EG(i, i) = £ + £v 2 EG(z, i)EU G (i) + ^(i) + 2 (z) - £ ]T e jp , (3.24) 

bl<« 

with 

\p\<n 

where 0i, 02 and e, p given by relations (3.18), (3.19) and (3.16). Let us denote 
the average EG(i, i) by g(i) and rewrite (3.24) in the following form : 

g(i) = £ + & 2 g(i)U g (i) + ^(i) + - £ Yl e ^ ( 3 ' 25 ) 

|p|<n 

Let us consider the solution {r(i), \i\ < n} of equation 

r(i) = £ + £v 2 r(i)U r (i), \i\ < n. (3.26) 

Given z G A,,, one can prove that the system of equation (3.26) is uniquely 
solvable in the set of iV-dimensional vectors {~F*} such that 

2 

ll^ll = sup\i\ <n \r(i)\ < (3.27) 

\lm{z)\ 

(see Lemma 6.1 of [ID])- Certainly, r(i) and g(i) (3.25) depend on particular 
values of z, n and b, so we shall use the notations r(i) = r n b (z; z) and g(i) = 
g nj h{i\ z). The following statements concern the differences : 

D n ,b(i; z) = g n , b (i] z) - r n>b (i; z), d n>b (i] z) = r n>6 (z; z) - w(z), 
where w(z) verifies equation (2.9). 
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Lemma 3.2 Given e > 0, there exists a positive integer L = L(e) such 
that for all sufficiently large b and n satisfying (2.3), inequality 

sup \d n>b (i; z)\ < e, z G A,,, (3.28) 

ieB L 

holds, with Bl given by (3.21). 

Lemma 3.3 If z G and (2.3) holds, then 

sup \D n , b (i,z)\ = o(l), n,6 — > oo, (3.29) 

|p|<n 



Lemma 3.1 follows from (3.28) and (3.29). Under the same conditions of 
Lemma 3.1, one can find L > L such that 

-PjT3^<I)--MI^ (3-30) 

Relation (3.30) follows from (3.22) and (3.27), a priori estimate 

sup|<7(i)| < — !-r, (3.31) 
\i\<n \lmz\ 

and the observation that l! has to satisfy condition ip(L — l!) < e. 

Proof of Lemma 3.2. Let us consider function Wi(z) = w(z) satisfying 
(2.9) that we rewrite in the following form similar to (3.26) : 

1 n 

Wi{z) = £ + £v 2 Wi(z)-^2b6 it Wt(z), \i\ < n. 

t=i 

Subtracting this equality from (3.26) , we see that d(i) = d ntb (i; z) is given 
by relation 

d(i) = tv 2 d(i)U r (i) + & 2 w{z)U d {i) + £v 2 w 2 (z)[P b + T{i)\, (3.32) 

where 
and 

Tn, b (i) = T(i) = \J2 ^(^) - I E K y (3 ' 34) 

\t\<n ' t& 
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It is clear that \P b \ = o(l) as 6 — > oo. Indeed, one can choose an even step-like 
function ipdif) , t G M such that 

fceN 

Then ^(t) < VK*)) V'dCO — * VKO as 6 ^ oo and the Beppo-Levy Theorem 
implies the convergence of the corresponding integrals of (3.33). 
Using equality 

we obtain relation 

= v 2 wr{i)U d {i) + v 2 w 2 r(t)[P b + T n>6 (i)], |z| < n, (3.35) 

where we write u> = w(z). This relation, together with estimates (3.27) and 
1^(^)1 < \Imz\~ 1 , implies that \d(i)\ < 2\Imz\~ 1 < 2. 

Given e > 0, let us find a number Q that 2 tp(t)dt < e. Setting r : = 
v 2 77~ 2 < 1/4, we derive from (3.35) such 

2 

sup < t[ sup + sup \T njb (i)\ +P b + 2e+ -], 

where we used the condition (2.1) and the estimate 

1 i Q — 7 /"°° 9 

j E \d(s)m s -^)<d mdt+- b 

s:\s-i\>Qb " ^ Q 

that follows from the monotonicity of ip(t). Now, it is easy to conclude that 

sup K*)| < £ r J '[ sup |T„ if) (2)|+P 6 ] + r L / Q sup|^)|+4e+^. 
ieB L 0<i<L/Q l*l<« ft 

(3.36) 

Let us choose such M that r M < e, then 



0<j<L/Q ieBL ~iQ ieB L - MQ 



^ r 7 sup |T n>6 (i)| < sup |T n>6 (i)| +2e. 
Finally, observe that 

2 ^ t /"°° 1 

sup |T n , b (i)| < - V ^(-) < 2 / V>(s)rf S + -. 

ieB L . MQ b t=n-(L-MQ)b ° J n/b-L+MQ 
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This inequality shows that (3.28) holds for sufficiently large L and 1 << 

b « n. o 



Proof of Lemma 3.3. Subtracting (3.26) from (3.25), we obtain relation 
for D(i) = D n>b (i,z), 

D{i) = iv 2 g(i)U D (i) + tv 2 D{i)U r {i) + ^(i) + 4> 2 {i) - i € *p- 

\p\<n 

This relation can be written as 

D = [I- W^Y 1 *), 
where we denote by W^- 9 ^ a linear operator acting on the vectors e = 

( e s)— n<s<n 

[W M e](i) = v 2 g(i)r{i) e(s)u(i,s) 

\s\<n 

and 

©o = ^{^i(o+02(o-^E^}- 

P =i 

It is easy to see that if z G A,,, then the estimates (3.27) and (3.31) imply 



\\W {9 ' r) \\ < j-^-rr < l- (3.37) 

11 " ~ \Imz\ 2 2 v ; 



sup ^{Uiii) + <h(i) -Z = z e ^ ( 3 - 38 ) 



Let us accept for the moment that 
^ 1 

in the limit n, b — > oo. Then, Lemma 3.3 follows from (3.37) and estimate 
(3.38). 

Now, let us prove (3.38). Using inequality (3.16), we obtain the first esti- 
mate concerning e« p , 

|V u r \ < £ ^ = 0(4=), (3.39) 

where we used condition (2.1) and 



|p|<n tea Jo 
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The term in (3.18) is estimated with the help of the elementary in- 
equality 

£|G(^)| 2 = ||Gel|| 2 < — (3.40) 

\p\<n 

which follows from (3.7) and the inequality ||G 2 (,2)|| < \Imz\~ 2 . Then, (3.40) 
implies 

|p|<n 

Using relations (2.1), (3.7) and identity 

E{^ } = E{/V}, 

we obtain that 

I^WI HyE V{G(i,i)G°(p,p)}4>C-^)\ < -^-E|C°(i,i)C/°(z)|. 

\p\<n 

(3.42) 

Thus, to prove (3.29), it is sufficient to prove inequalities (3.39), (3.41), 
(3.42) and the estimate 

supE|G°(M)C/£(i)| =o(l), (3.43) 

|i|<n 

We prove estimate (3.43) in section 4. Assuming that this is done, one can 
say that Lemma 3.1 and relation (2.11) are proved. 



4 Proof of the main Theorems 

In this section, we study the variance Var{g nt i ) (z)} and complete the proof 
of Theorem 2.3. Finally, we prove Theorem 2.2 and Theorem 2.1. 

4.1 The variance and the proof of Theorem 2.1 

In section 3, we proved relation (2.11) assuming that the asymptotic re- 
lation (3.43) holds. To complete the proof of Theorem 2.3, we follow mainly 
two steps. The first step is to remark that the estimate (3.43) and (2.12) are 
consequences of the fact that the variance of the diagonal entries G(i,i;z) 

Var{G(i,i;z)} = E\G°(i, i; z)\ 2 
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vanishes as n, b — > oo. The second step is to prove the following two relations 
that concern the moments of diagonal elements of G. 

Lemma 4.1. If z G A v (2.10) and fa < oo (2.4), then the estimates 

supE|f/°(2)| 2 = o(6- 3/2 ) (4.1) 

|i|<n 

anc? 

supE|G (i,i;z)| 2 = o(fe- 1 / 2 ), (4.2) 

|j|<n 

are £nte in the limit n,b — > oo. 

It it easy to show that estimate (3.43) follows from (4.2). Finally, (2.12) 
follows from (4.2) and inequalities 

Var{^)} = E|^)-E(^ t (^))| 2 <E{^ £ |G°(p,p)G°( a , s )|} 

|s|,|p|<n 

<^ E V^PMv^PMP) < sup E(|G°(p,p)| 2 ). (4.3) 

M,|p|<n ^ n 

Theorem 2.3 is proved, o 

We close this subsection with the proof of Lemma 4.1. 

Proof of Lemma 4-1- Let us denote G± — (H — ^i) _1 , G 2 — (H — z 2 )~ 1 
with Zj G A v , j = 1,2. Regarding the average E{G5(i, i)G 2 (p, p)}, we apply 
to G 2 the resolvent identity and we use (3.2), to get 

E{G?(M)G 2 (P,P)} = -6 E E{GS(i,0G 2 (p,a)}. 

\s\<n 

For each pair expression i)G 2 (p, s) represents a smooth function 

of H(s,p) and its derivative is bounded because of (3.6) and (3.7). In particu- 
lar \D 2 sp [G\{i,i)G 2 {p, s)}\ < ci{\Imzi\~ l + \Imz 2 \~ l ) A where c\ is a constant 
and Dip = d r /dH(s,p) r . 

According to (2.2), (1.5) and condition fa < oo (2.4), the third absolute 
moment oiH(s,p) is of the order b~ 3 ^ 2 . Thus applying (3.4) to E{Gl(i, i)G 2 (p, s)} 
with q — 1, one obtains 



EG°^^)G 2 (p,p) = Z 2 v 2 E{G ^^)G 2 (p,p)U G M}H2V 2 £ E{G?(*, i)G 2 (p, s) 2 «(a,p)} 

|s|<ra 
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+2&v 2 s)Gi(p,i)G 2 (p, s)u(p, s)} - ^E{G°(M)G 2 (p,p) 2 

\s\<n 

+G 1 (i,p) 2 G 2 (p,p)}i;(0)-^J2 T ^ 

\s\<n 

where 

^^{suppy^^ (4.4) 
with a constant c, 

^(5=2) x ^ 
u(s,p) = — - — and U G {p) = } j G(s, s)u(s,p). 

\s\<n 



Let us introduce variables 

&V 2 U E (G 2 )(P) 



qM = i-eJL^toV M ~ (45) 



Using (4.5) and identity 

E{fg} = E{fg°} + E{f}E{g}, (4.6) 

we get equality 

W%(i,i)G 2 (p,p) = q 2 (p)v 2 E{G° l (t,t)G 2 (p,p)U G2 (p)} 

+q2(p)v 2 nG° 1 ^t)G 2 (p,s) 2 u(s,p)} 

\s\<n 

+2q 2 (p)v 2 E{G 1 (i,s)G 1 (p,i)G 2 (p,s)u(p,s)} 

\s\<n 

-^^E{G\{i,i)G 2 {p,p) 2 + G 1 (t,p) 2 G 2 (p,p)}m - q 2 (p) J2 r sp - (4.7) 

\s\<n 

Multiplying (4.7) by u(i,t) and summing over i, we get relation 
EU° Gl (t)G 2 (p,p) = q 2 (p)v 2 E{U Gl (t)G 2 (p,p)U G2 (p)} 

+q 2 (p)v 2 V{U Gl {t)G 2 (p,s) 2 u{s,p)} 

\s\<n 

+2q 2 (p)v 2 E{Gi(i, s)G 1 (p, i)G 2 (p, s)u(p, s)u(i, t)} 
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\i\<n 

~^<l2{p)Tspu{i,t). (4.8) 

i,s 

Regarding Gi(p, .)«(., t) and G 2 (., s)u(s, .) as n— dimensional vectors, one 
gets 

| J2 E ( G i(^ s)Gi(p, i)G 2 (p, s)}«(p, s)«(i, t) | 



< \\Gi\\ E |Gi(p, <)«(<, 0l 2 /E |G 2 (p,s)«(s,p)| 2 . (4.9) 

Inequalities (3.7) and (3.40) imply that the right-hand side of (4.9) is bounded 
by b~ 2 r]~ 3 . Regarding the last term of (4.8) and using (4.4) and (3.7), we get 



i E ^,t)i<-^/i3[E^^)][E M ^ s )] = o(4f)- (4.io) 

|i|,M<» 16V Vb \i\<n \s\<n V b 

The first term of the right-hand side of (4.8) can be estimated by 



\q 2 (pyE{U Gl (t)G 2 (p,p)U° G2 (p)}\ < -^—- 2 /sup EI^COI 2 . /sup E|t/g 2 (i)|*, 

\imz 2 \ z Y |j|< n y|j|<n 

(4.11) 

where we used estimates (3.7) and 



MP)\ < (4.12) 



Also, we use (3.40) and (4.12) to see that 



\ q2 ( P y e E{f/° 1 (t)G 2 (p, S ) 2 M ( S , P )}i < — ^—- 3 isu P {E\u Gi m 

\s\<n 1 2| V |l| " n 

(4.13) 



|s|<ra 

and 



|^^E{^ i (t)G 2 (p,p) 2 }^(0)+^^E{G 2 (p,p) E GiCvp) 2 ^*, 0}^(0)| 

|i|<n 

* ^l^™^ + ^^W - <414) 
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Now, multiplying (4.8) by u(p, r) and summing over p, and taking G\ = 
Cx2 and i = r, we obtain a relation that together with (4.9), (4.10), (4.11), 
(4.13) and (4.14) implies the following estimate for the variable M 12 = 
sup mn {E\U Gl (i)\*} : 



,2 



M 12 < ,/ , 2 M 12 + 4 1 V / ^+ 4|, 

with and A 2 are constants. Regarding this inequality, it is easy to show 
that M12 = 0(b~ 3/2 ). This proves (4.1). 

In order to prove (4.2), we go back to relation (4.7) with G\ = G 2 , p = i. 
Applying (3.7) and (4.12), we obtain estimates 



V 2 AZTTTT 7~7~" t> 2 



sup |G?(M)| 2 < -W^I^OI 2 ) + l^iC*. 

|i|<n 7 7 u 



S 



|2 



\s\<n 

2v 2 



' |s|<n ' ' |s|<ra 

with 7] = 2v + 1. Regarding the last term of (4.15) and (4.4), we obtain 
inequality 

|s|<n |s|<n 

where c 2 is a constant. Applying (3.40), (4.1) and (4.16) to inequality (4.15), 
we obtain (4.2) and we are done, o 

4.2 Proof of Theorem 2.2 and Theorem 2.1 

Proof of Theorem 2.2. It is easy to show that relation (2.11) together with 
(2.12) implies convergence in probability : 

P - lim g n (z) = w(z), z G A„. 

n^oo 

By definition (2.6), we rewrite this convergence in the form 

{ / ^f^}^ -A { / d A z e A,. (4.17) 

Now we will prove that (4.17) is true for all non real z. Representation (4.17) 
for g n (z) and w(z) implies that the function z i— > g n {z) — w(z) is analytic 
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and uniformly bounded on any compact set T C A,,. Thus, given e > 0, there 
exists a finite set z lt Zk G T such that 

max \g n (z) - w(z)\ < e + max \g n (zj) - w(z j )\. 
zer j=i,..,k 

Let {ni} a subsequence from {n}. We have g ni (zi) — * w \ z x)i then there 
exists a subsequence {n x } C {n\\ such that g n '{z\) -^4 w(zi). Therefore, 

the sequence {%} contains a subsequence {n 2 } where g n ' 2 (z 2 ) 117(22) and 

g n ' 2 (zi) Now we see that there exist a subsequence n = n' k C ... C 

n l common to all 21, such that 

max. k»'(*i) - ^ °- 

j— i,..,/e 

As a result, we get consequence 

lim {max |<7 n '(2) — 7/7(2)!} = a.s.. 

n — >+oo ^GT 

Now we obtain g n >{z) — ^* 7/7(2) with {n } is common for all 2 G T. Since 
g n >{z) and 7/7(2;) are analytic functions in C \ R, and g n >(z) 7/7(2) for 
zeT, then this convergence is true for all 2 6 C \ R. Thus, we rewrite (4.16) 
as 

d(T n',b'(^ H n') , a.s. , /" tfo"(A) 



The sequence {ni} contains a subsequence {n } such that (4.18) holds, 
which implies convergence in probability (2.9) (see Lemma in paper |11|). o 

Proof of Theorem 2.1. It is well known that the set of linear combinations 
of {xrj"j z j ^0} is dense in C (R) (set of counting functions), then (4.18) 
implies that 

{ I ^(Xjda^iX,^)}^ ^ { [ V (X)da(X)} (4.19) 
Jn Jr 

for all ip G C (R). 

The sequence {rii} contains a subsequence as {n } such that the conver- 
gence (4.19) is true. We use the same arguments as in the proof of Theorem 
2.2, the convergence in probability (2.5) gives and finishes the proof of Theo- 
rem 2.1. o 
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5 Estimates of the Variance 



In this section, we obtain more estimates of the variance under more 
restrictive conditions on the moments of a(i,p) and the function ip(t) given 
by (2.1). Finally, we provide a proof of the cumulant expansions formula. 

5.1 First Estimate 

Theorem 5.1. Assume that E{a 2l+1 (i,p)} = with I = 0,1,2 and 

E{a(i,p) 2 } = v 2 , E{a(i,p) 4 } = 3v 4 , E{a(i,pf} = 15v 6 , \p\ < n, 

fij < oo (2.4), then the estimate 

Var{g n , b (z)} = 0{\) (5.1) 

holds for large enough n and b and for all z e A v . 

Proof of Theorem 5.1. We prove Theorem 5.1 by using the following esti- 
mates of the moments of the diagonal elements of the resolvent G 

Lemma 5.1. Under conditions of Theorem 5.1, the estimates 

su P E(|f/°(^)| 2 ) = 0(i) (5.2) 

|i|<7» 

and 

su P E(|G(w)| 2 )=0(i) (5.3) 

\i\<n 

hold in the limit n, b — ► oo and for all z <E A v . 

Now it is easy to show that Theorem 5.1 follows from inequality (4.3) and 
estimate (5.3). 

Proof of Lemma 5.1. We start with (5.2). Let us denote G\ = (H — zi)^ 1 , 
G 2 = (H—z 2 )~ 1 with Zj G A v , j — 1, 2. Consider the average E?7^ l (t)G ? 2 (p,p) 
and use (3.2) to obtain 

n 

EU° Gl (t)G 2 (p,p) = -Cj2 E { U k(t)G2(p,s)}. 

s=l 
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For each pair (s,p) U G (t)G 2 (p, s) is a smooth function of H(s,p) and its 
derivative is bounded because of (3.7). In particular \D^ [U Gl (t)G2(p, s)]\ < 
cz{\Imz\\~ l + \Imz2\~ 1 ) 8 for some constant C3 and D r sp = d r /dH r (s,p). 

According to (2.2), (1.5) and condition £l 7 < 00 (2.4), the 7-th absolute 
moment of H(s,p) is of order b~ 7 ^ 2 . Thus, applying (3.4) to EU G (t)G2(p, s) 
with q — 5, one obtains 

E{U Gl (t)G 2 (p,p)} = -6 E K ^{DlM^)G2(p,s)]} 

\s\<n 

E y E {^[^(^2(f^)]} -6 E ^{^[^(^(p,*)]} 

|s|<n |s|<n 

-« 2 E 

|s|<n 

where 

|e«| < csupl^y^^^^^lir^Eda^DE^). (5.4) 

Regarding the terms involving the first derivative and using (3.6), the defi- 
nition of K 2 (3.13), we get 

E{U° Gl (t)G 2 (p,p)} = Z2V 2 V{U Gl (t)G 2 (p,p)U G2 (p)} 
H2V 2 J2v{U Gl (t)G2(p,s) 2 }u(p,s) 

\s\<n 

+2& 2 E E{G l (t, s)G l (p, i)G 2 (p, s)}u(p, s)u(i, t) + R(p), 

s,i 

where 

Rip) = -^fE{U Gl (t)G 2 (p,p) 2 + E Gi(i,p) 2 G 2 (p,p)«(<,*)}^(0) 

\i\<n 

|s|<n |s|<n 

E £>■ 

|s|<n 

This term does not contribute to the estimate (5.1). 

Applying (4.6) to the first term and using the definition of (72 (p) (4.5), 
one gets 
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E{U° Gl (t)G 2 (p,p)} = q 2 (p)v 2 E{U Gl (t)G 2 (p,p)U° G2 (p)} 
+q 2 (p)v 2 V{U Gl {t)G 2 (p,s) 2 u{s,P)} 

\s\<n 

+2q 2 (p)v 2 E{Gi(<, s)G l (p, i)G 2 (p, s)u(p, s)u(i, t)} + ^-R(p). (5.5) 

st i2 

Inequality (4.9) implies that the third term of (5.5) is of order 0{b~ 2 ) as 
n, b — > oo. 

Let us estimate each part of (5.5) with the help of U G . Using (3.40) and 
inequality (4.12), we can estimate the first term of the right-hand side of 
(5.5) by 



v 2 



\q 2 (p)v 2 E{U Gl (t)G 2 (p,p)U G2 (p)}\ < — — 2 supE\U Gi (i)\isupE\U G2 (i)\i. 

\lmz 2 \ z Y |j|< n |j|<„ 

(5.6) 

Also we can write that 

\q 2 (p)v 2 £ E{U Gi (t)G 2 (p,s) 2 u(s,p)}\ < — ^— /supjEl^WP}. 

| s |<n b \ LmZ2 \ V 1 ^" 

(5.7) 

If one assumes for a while that 

sup \R( P )\=od + W° 1 (*)|)}, (5.8) 

\p\<n 

holds, then (5.2) will follow. Indeed, multiplying (5.5) by u(p, r) and summing 
over p, and taking G± = G 2 and i = r, we obtain a relation that together 
with (4.9), (5.6), (5.7) and (5.8) imply the following estimate for variable 

M 12 = sup mn \U Gi (t)\ 2 : 

M12 < ^M 12 + -±y/Afa + 

with A 3 and are constants. Regarding this inequality, it is easy to show 
that M12 = 0(6~ 2 ). This implies (5.2). 

Let us prove (5.8). Using inequality (5.4), we get 

\Y,^\<j^Y,^p) = °(w^ (5 - 9) 

\s\<n V /; \s\<n 
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for some constants c and c 3 . By the definition of the cumulants (see subsection 
3.1), we obtain 



K 4 (Y) = E(F 4 ) - 3E(F 2 ) 2 = -A sp , (5.10) 
where Y = H(s,p), A ip = 3^(^){1 - V(^)} and 

K 6 (Y) = E(F 6 ) - 15E(F 4 )E(F 2 ) - 10E(F 3 ) 2 + 30E(F 2 ) 3 = ^Q sp (5.11) 

with Q ip = 15^){1 - ff^) + §^) 2 }- 

Using (3.6) and (3.7), we obtain that sup Mp \D r sp [U^ l (t)G 2 (p, s)}\ = 0(1) 
with r = 1,2,... Which gives 

|s|<n |s|<n 

where Mi is a constant. 

One can estimate the terms with the third derivative by using similar 
estimates in relations (4.9) and (5.6) : 

^ E ^{DU U k(t)G2( P , *)]} = 0{1 + - b E(\U° Gl (t)\)}- (5-13) 

|s|<ra 2 

Indeed, E{D^ p [UQ i (t)G 2 (p, s)}} contains 44 terms of the form 

^2 E i G n («i> Pi)G r2 (a 2 , I3 2 )G 13 (a 3 , /3 3 )C 74 (a 4 , /3 4 )G 75 (a 5 , /3 5 )}u(i, *)» 

|j|<n 

and 3 terms of the form 

72(^2, P2)G ri («3, /?3)G 74 (o;4, /34)G 75 (q;5, /3s)}li(i, t), 

|i|<n 

where 7^ = 1 or 2 and otj G {i,s,p} and /3j G {i, s,p}. These terms can be 
gathered into three groups. In each group the terms are estimated by the 
same values with the help of the same computations. 

We give estimates for the typical cases. Using (2.1) and (3.40), we get for 
the first terms of the first group : 

\J2MH(s,p))E{G 1 (i,p) 2 G 1 (s,s)G 1 (p,s)G 2 (p,s)}u(i,t)\ 

i,s 
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q4 o..4 
|i|<n |s|<n 

For the terms of the second group, we obtain estimates 
\^ s K 4 E{G 1 (i,s)G 1 (i,p)G 1 (p,p)G 1 (s,s)G 2 (p,s)}u(i,t)\ 

< p E E {l G i(^ ^0°- s ) IMP, S M*> *) 



^ 6i^i3,/E |Gi(p,iMt,t)| 2 I^MM^I 2 = O(i). (5.15) 

Finally, for the terms of the third group, we get inequalities 

| ]T K 4 E{G 1 (i, i)°G 2 (p, s) 2 G 2 (p, p)G 2 (s, s)}u(i, t) \ 



l,S 

,4 



Q 7 ,4 1 

-^ E(|f/ ^ (t)l){ S M(s ' p)} = 0( 6 E(|f/ ^ (t)l)) - (5 ' 16) 

|s|<n 

Gathering all the estimates of 47 terms, we obtain (5.13). Finally, notice 
that the term 

I^El^^^fe^+^d^p) 2 ^^,^^^)}^^)! 

\i\<n 

is estimated using (4.14). Then (5.8) follows from (4.14), (5.9), (5.12) and 
(5.13). o 



Let us prove (5.3). Using (3.2), one has 

EG?(M)G 2 (P,P) = £ E{GS(M)G 2 (p,a)}. 

|s|<n 

For each pair (s,p) Gl(i,i)G 2 (p, s) is a smooth function of H(s,p) and its de- 
rivative are bounded because of (3.6) and (3.7). In particular \D* [Gi(i, i)G 2 (p, s)]| < 
c^dlmzil" 1 + llmz^ 1 ) 6 with a constant c 4 . 

According to (2.2), (1.5) and condition /t 5 < 00 (2.4), the 5-th absolute 
moment of H(s,p) is of order 5~ 3/2 . Thus applying (3.4) to E{G°(i, i)G 2 (p, s)} 
with g = 3, one obtains 

EGl(i,i)G 2 (p,p) = ^v 2 E{G° 1 ^t)G 2 (p,p)U G2 (p)} 
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H2V 2 E{G° 1 ( i , i )G 2 (p,s) 2 u(s,p)}+^2V 2 E{G 1 ( i ,s)G 1 (p, i )G 2 (p,s)u(p,s)} 

\s\<n |s|<n 



6 

|s|<n M<n 



E fv{D%[G^)G 2 (p,s)]} -6 £ eg) 



with 



|e£?| < C {sup|^ 4 p [G?(^,.)G 2 (p, S )]|}r 5 / 2 E(|a ps | 5 )E(^ 



Using the identity (4.6) and the definition of q 2 (4.5), we write 
EG?(M)G 2 (P,P) = g 2 (p)t; 2 E{G?(^^)G 2 (p,p)f/° 2 (p)} 

+g 2 (i^ 2 £ e{g°(m)g 2 (p,s) 2 «(s,p)} 

|s|<n 

+2g 2 (p)t; 2 £ E{Gi(i,s)Gi(p,i)G 2 (p,s)«(p,s)} 

|s|<n 

^ i/Mt;2 E{G;^^)G 2 b,p) 2 + G 1 (t,p) 2 G 2 (p,p)^(0) 



b 

-<fe(p) E ^E{Dj,[GS(i,i)G2(p,s)]}-(b(p) E ( 5 - 17 ) 
It is easy to show that 



|s|<ra 



with M 2 some constant. Inequality (3.7) implies that Dl p [Gl(i, i)G 2 (p, s)] = 
0(1) with r = 1,2, .., whence one obtains inequality 



| K 4 E{D%[G l ( i , i )G 2 (p,s)]}\ < M ( s >^' ( 5 - 19 ) 

|s|<n |s|<n 

for some constant M 3 . Finally, it is easy to show that (3.7) implies estimate 

\ q -^^E{G^)G 2 (p, P r + G 1 (i,p) 2 G 2 (p,p)}^(0)| = 0{\). (5.20) 

Now, (5.3) follows from relation (5.17), inequality (3.40) and estimates 
(5.2), (5.18), (5.19) and (5.20). The estimate (5.3) is proved, then Lemma 
5.1 follows and finishes the proof of Theorem 5.1. o 
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5.2 Second Estimates for the Variance 

In this subsection, we derive estimates of the variance under more strong 
conditions on H n and ip. 

Theorem 5.2. Suppose that E{a 2l+1 (i,p)} = with I = 0, 1,2, 
E{a(i,p) 2 } = v 2 , E{a(i,p)*} = 3v\ E{a(i,p) 6 } = 15w 6 , \p\ < n 
and 

/tio = sup E\a(i,p)\ 10 < oo and / a/ ip(t)dt < oo, 



t/ien £/ie estimates 



\i\<n 

and 



sup |E«£(z)G°(M)l = + jyVar^)) (5.21) 



VarK(z)} = O(p) (5.22) 



/zoic? /or /an/e enough n and b and for all z <E A v . 

Proof of Theorem 5.2. We start with (5.21). The proof follows the lines 
of the proof of (5.3). Regarding (5.17) and summing it over i and using the 
notation g n (zi) — 9i, we obtain that 

NV{g°G 2 (p,p)} = ^g 2 (p)u 2 E{^G 2 (p,p)C^ a (p)} 
+Nq 2 (p)v 2 EKG 2 (p, s) 2 u(s,p)}+2q 2 (p)v 2 £ E^i, s)G 1 (p, i)G 2 (p, s)u(p, s)} 

\s\<n s,i 

_^! E {iv^G 2 (p,p) 2 + G 1 (t,p) 2 G 2 (p, P )}m 

\i\<n 

- n ^(p) e ^{u^G^^D-fe^x;^- ( 5 - 23 ) 

|s|<n s,i 

Using (5.2), (3.7), (3.40) and (4.12), it is easy to show that 

AT , 

Nq 2 (p)v 2 E{glG 2 (p,p)U« G2 (p)} = 0(- v / Var^)), (5.24) 
Nq 2 (p)v 2 £ EKG 2 (p, S ) 2 M ( S ,p)} = O(yv^Var^O), (5.25) 

\s\<n 
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< 



|2g 2 GV E i G ^ S ) G ^ l ) G ^ S ) M G°> s ))l 

s,i 



\i\<n y |s|<ra 



and 



^^-E{NglG 2 {p,pf + G 1 ( l , P ) 2 G 2 (p,p)M0) 



\i\<n 



= 0(^v/\M^) + i). (5.27) 

Now, regarding Dl p [glG 2 (p, s)] with r = 1, 2, .. and (5.31) (see Lemma 5.2 at 
the end of this subsection), it is easy to see that 



s,p 

Then we obtain 



S np\D: p [g 1 G 2 (p,s)]\=O(^ + \g° 1 \), r = l,2,... (5.28) 



\Nq 2 (p) K^{DlMG 2 {p,s)]}\ <^Y1 ^F\V{D 3 s MG2(p,s)]}\ 

\s\<n ' \s\<n 

,N r 1 



= 0( T [- + v / Var^y]). (5.29) 
Now we give more details for the remainder 

|eg)| < c{K 4 E\H(s,p)D 4 sp [G° 1 (i,i)G 2 (p,s)}^\ + K 2 E\H 3 (s,p)Dt p [G° 1 (i,i)G 2 (p, s)}^\ 

+E\H 5 (s,p)D%[G°^^)G 2 (p,s)}^\}, (5.30a) 
where for each pair (p, s) 

[G° l (z,t)G 2 (p,s)}^ = {G^}° sp (i,i;z 1 )G^(p,s;z 2 ), v = 1,2,3 
and G%\ Zj ) = {H%> — Zj) 1 , j — 1, 2 is the resolvent of the real symmetric 



matrix : 



Tr(u) (r -\ _ f H(r,i) if (r,i) ^ (s,p) 



y) (s,p) if (r,i) = (s,p) 
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with \H^(s,p)\ < \H(s,p)\, v = 1,2,3, (see Lemma 5.3 at the end of this 
subsection). Regarding (5.30a) and summing it over % and s, we obtain 



E l e 2 } l ^ ^ c ^^{■K4E\H(s,p)D / j ip [giG2(p, s)]^\ + K2E\H 3 (s,p)D^ p [g^G 2 (p, s)]^\ 

s,i \s\<n 

+E\H 5 (s,p)D%[g 1 G 2 (p,s)f'>\} (5.306) 

Hereafter, we use the notation : for each pair {p, s) and v — 1,2, 3, let 
= H be the matrix defined by 



H(r,i) 



H(r,i) if (r,i)^(s,p) 
H(s,p) if (r,i) = (s,p) 



with \H(s,p)\ < \H(s,p)\, the resolvent G^p(zj) = Gj and its normalized 
trace by gj, j = 1, 2. To derive estimates of e^ 2 \ we need the following Lemma. 

Lemma 5.2. If z G A,, (2.10) and for large enough n andb, the estimates 
D r s M(z)}=0(^), r = l,2,.. (5.31) 

and 

Var(g n (z)) = 0{Var(g n (z)) + ^}, (5.32) 
are true in the limit N, 6 — > oo. 

We prove Lemma 5.2 at the end of this section. 

Keeping in mind the expression (5.306) and using (5.28), (5.31) and (5.32), 
we obtain inequalities 

E m\s, P )DtM^ s )}\ < ^ E ^ + p E n\H{s,P)m\] 

\s\<n N<n \s\<n 

|s|<ra |s|<n 

^ ^ E — — + -pT^V Var{ 5l } ^ — ^- 

|s|<ra 

1 r 1 



= °{^[^ + V / VarTM]} ) (5-33) 
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J2K,E\H( s ,p)Dt p {g°G^s)}\ 

\s\<n 

- ]V6 3 / 2 ^ b b 3 / 2 ^ b 

|s|<n H<n 
1 r 1 



and 



= 0{^[- + v / V^M]} (5.34) 
5]^ 2 E|^ s ( a ,p)£>J,{^S(? a )}| 

\s\<n 

3^3 v Vgg) 3^ v /^y/Va7{M ^ Vgg) 

|s|<n H<n 

1 r 1 



= 0{^[- + v / Va^]} (5.35) 

with some constant f3 and fi r = sup ip E|a(i,p)| r < oo, r = 1, .., 10. Then 
(5.33), (5.34) and (5.35) imply the following estimate 

\J2 e( S\ = o{N[^ + Vy^i))]} 

s,i 

= O[JV(^ + ^0MS))]- (5 - 36) 

Returning to (5.23), using definition N = 2n + 1 > n and gathering 
estimates (5.24), (5.25), (5.26), (5.27), (5.29) and (5.36), we get (5.21). 



Now, let us prove (5.22). Regarding (5.23) and summing it over p with 
G 2 = Gi = G and using (5.21), we obtain estimate 

Var{g(z)} < A{± + ±y/Var(g(z))). (5.37) 

where A is a constant. Regarding this inequality, it is easy to show that 
Var{g n (z)} = 0(b~ 2 ). Relation (5.22) is then proved and ends the proof of 
Theorem 5.2. o 
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Proof of Lemma 5.2. We start with (5.31). Consider the normalized trace 
9n{z) = g{z) = ^Y.\k\<n G ( k ' k )^ witn z e \- Regarding (3.6), we obtain 

Dl p {g°(z)} = 1 £ ^{G(k,k)} = -| E G(k,s)G(k,p) 

\k\<n \k\<n 



~G"(»,rt. 



Using (3.7), we get 



Dl p {g°(z)}\ < 



sp ^ v ,j, _ N \ Imz \2- 

Regarding the second derivative, one obtains 

D%{g\z)} = l{2G 2 (s,p)G(s,p) + G 2 (p,p)G(s,s) + G 2 (s,s)G(p,p)}, 

using (3.7), we obtain D^ p {g°(z)} = 0(-^). Then it is easy to show that 
Dr s P {g°(z)} = 0(jj), r = 1,2... Estimate (5.31) is proved. 



Now let us prove (5.32). Using the resolvent identity (3.1) for two hermi- 
tian matrices H and H , we obtain 

g = Uz) = ^Yl k ) = E G ^ V+jj E r ) fi-H\ (r, i)G(i, k) 

\k\<n \k\<n k,r,i 

= g+±Tr(GG5 H ) (5.38) 



if (r,i)^(s,p), 

H(s,p) - H(s,p) if (r,i) = (s,p). 



with 

[6 H ](r,i) = [H-mr,i) = 
Then 

Tr(GG5 H ) = Y,[GG](r,i)[5 H ](r,i) = [GG](s,p)[H(s,p) - H(s,p)\. (5.39) 

r,i 

Let us return to relation (5.38). It is easy to see that 

Var(^) < 2Var((?) + ^Var(Tr{GG5 H }). (5.40) 
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Then assuming that 

Var(Tr{GG6 H }) = 0( J) (5.41) 

holds, (5.32) follows from (5.40) and (5.41). Let us prove (5.41). It is clear 
that \H sp — H sp \ < 2\H sp \. Using (3.7) and (5.39), we obtain 

Var(Tr{G«5fl}) < jj^E(\H(s,p)\ + E(\H(s,p)\)) 2 

< tylt^ miaspl){1 + ^( £ ^))lV( £ i^) + E(|a sp |) 2 ^(^) 2 (l - ^( S — 

= 0(i). o 
Finally, we prove (5.30) for the remainder e^. 

Lemma 5.3. (The cumulant expansions method with q=3). Let us 

consider the family {Xj, j = 1, ..,m} of independent real random variables 
determined on the same probability space such that E{|X,-| 5 } < oo, j = 
1, .., m and 

E(Xj) = E(Xf) = 0, j = l,..,m. 

If F(ti, .., t m ) is a complex-valued function acting on vectors (with compo- 
nents as a real variables) in a m— dimensional space such that its first 4 
derivatives are continuous and bounded, then, for all j , we obtain (3.4) with 
q = 3 and the remainder 6% is given by 

es = -f EWf (Y 2 )} - §E{X 3 /jVi)} 

+ ^E{X 5 /f(>o)}, (5.42) 

where \Y V \ < \Xj\, v = 0,1,2, K r is the r-th cumulant of Xj and fj is a 
complex-valued function of one real variable such that 

fj(Y) = F(Xi, ..,Xj-i, Y, Xj+i, ..,X n ) 

and f^ is its 4—th derivative. 

The cumulants can be expressed in terms of the moments. If [i r = E(XJ) 
with j = l,..,m, one obtains : K\ — [i\ — 0, K 2 = [i 2 , K 3 = 0, K 4 = 
Hi -3(4. 
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Proof of Lemma 5. 3. To simplify the notation we write X = Xj for a real 
random variable such that E{|X| 5 } < oo and E(X) = E(X 3 ) = 0, and we 
denote / = fj. 

Now, consider a function / of the class C 4 and apply Taylor's formula to 
/, to obtain relation 



3 ~\ t~ f v4 



= E 7r/ (r) (o) + ^r/ (4) (n), \y \ < \x\. 

r=0 

Multiplying this equation by X and taking the mathematical expectation of 
both sides, we obtain equality 

E{Xf(X)} = E(X 2 )E(/'(X)) + E(|^)/( 3 )(0) + E{^f^(Y )}. (5.43) 

Now we apply Taylor's formula to /', and write that 

/'(X) = /'(0) + x/( 2 )(o) + |^)(0) + ^fWfr), 

with | li| < \X\. Taking the mathematical expectation of both sides and 
multiplying by E(X 2 ), we obtain relation 

E(X 2 )E{/'(0)} = E(X 2 )E{/'(X)} - [E( ^ 2)]2 E{/( 3 >(0)} 

-E(X 2 )E{^/( 4 )(F 1 )}. (5.44) 
Applying Taylor's formula to f( 3 \ we see that 

E{/( 3 )(0)} = E{/( 3 )(X)} - E{Xf^\Y 2 )} 
with |l2 1 < Then we can write that 

E(X 4 )E{/®(0)} = E(X 4 )E{/^(X)} - E(X 4 )E{X/( 4 )(F 2 )} (5.45) 
and that 

E(X 2 ) 2 E{/( 3 )(0)} = E(X 2 ) 2 E{/( 3 )(X)} - E(X 2 ) 2 E{X/( 4 )(F 2 )}. (5.46) 
Now, substitute (5.46) in relation (5.44) and use (5.45), to rewrite (5.43) as 

B{Xf(X)} = E(X 2 )E(/'(X)) + |iE{/( 3 )(X)} 
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+ i. E {X 5 / (4) (^o)} - ff E W (4) ^)} " ff E W (4) (^i)} 

with K 4 = [E(X 4 ) - 3E(X 2 ) 2 ] and K 2 = E(X 2 ). This gives relation (5.42). 

o 
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